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Abstract. In this paper we study the late evolution of model galaxies after a single episode of star formation of different 
durations. The aim of the paper is to discover the timescale needed to refill with cold gas the center of the galaxy. This timescale 
strongly depends on the amount of gas initially present inside the galaxy and ranges between 125 and 600 Myr. A H i hole can 
therefore survive several hundred Myrs after the last SNII has exploded. If, as a consequence of the refill of the center of the 
galaxy, a second episode of star formation occurs, it pollutes the surrounding medium in a very short timescale (of the order of 
10-15 Myr), at variance with what happens if the center of the galaxy is still occupied by hot and tenuous gas. 
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1. Introduction 

After an episode of star formation (hereafter SF), the energy re- 
' leased by SN explosions and stellar winds produces large cav- 
ities of hot gas (visible in H i as holes), surrounded by a cold 
H I shell. However, in many dwarf irregular galaxies (hereafter 
, dirr) there are H i holes not associated with young star clusters 
■ (e.g. MlOl, Kamphuis et al. 1991; Holmberg II, Rhode et al. 
1999; LMC, Kim et al. 1999). 

The energy input rate after an episode of SF declines with 
time, therefore the hot cavity looses pressure and the shell 
tends to recede towards the center of the SF region. This re- 
. fiU process strongly depends on the considered energy sources. 
' Superbubbles solely produced by SNell experience a buoyancy 
of the hot gas at the end of the SF phase and the central re- 
gion can be completely refilled with cold gas. This process 
occurs in a timescale of the order of a few 10^ yr (D'Ercole 
& Brighenti 1999, hereafter DB99). Also Hunter & Gallagher 
(1990) speculated on the final fate of a supershell, based on 
simple dynamical arguments and estimated a refill timescale of 
the order of ~ 5 x 10^ yr The assumption of an energy source 
coming also from SNela (which should indeed be numerous in 
irregular galaxies; see e.g. Mannucci et al. 2005) changes the 
thermodynamical behaviour of the gas in the center of galax- 
ies. For instance, Recchi et al. (2001) found that in some cases 
the break-out of the galaxy occurred only as a consequence of 
the energy input from SNela, while the SNell, exploding in a 
colder and denser medium, could radiatively loose a very large 
fraction of the initial explosion energy. 

In some cases, the overpressurized regions produced by the 
ongoing SF develop galactic winds. Galactic winds are defined 
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as outwards-directed flows of gas, whose expansion velocities 
reach the escape velocity, therefore supposed to definitely leave 
the parent galaxy. In stratified media (i.e. if the density profile 
in the vertical direction is p, oc 2^'), a galactic wind develops 
for profiles steeper than x=2 (Koo & McKee 1992). However, 
the final fate of the outflowing gas depends on poorly con- 
strained details about the environment surrounding the galaxy 
and about interactions with other galaxies and it is hard to dis- 
tinguish, both theoretically and observationally, galactic winds 
from normal outflows. However, the refill of a superbubble is 
not affected at all by the final fate of the outflowing gas, there- 
fore this distinction is not important in our work. 

The refill of the center of the galaxy changes also drasti- 
cally the chemical evolution of dirrs. If previous episodes of SF 
carve a very large cavity, the release of metals from dying stars 
occurs in a very hot medium. Under these conditions, the cool- 
ing timescale of the newly produced metals is very large and, 
in the presence of a galactic wind, these metals can be directly 
carried out of the galaxy (Recchi et al. 2005; hereafter Paper I). 
Consequently, the chemical enrichment of the galaxy is not in- 
fluenced by the ongoing production of metals and the chemical 
composition only reflects the enrichment from the first episodes 
of SF. However, if the preceding generations of stars do not pro- 
vide enough energy to create a galactic wind and to produce a 
large cavity of hot gas, the impact of the following episode of 
SF can be significant (Recchi et al. 2002). The same occurs if 
the gap between two episodes of SF is large enough to allow 
the refill of the center of the galaxy with cold and dense gas. 

In this paper, we simulate the evolution of a galaxy after a 
single episode of SF of different duration, in order to establish 
what is the typical refill timescale when the energy input of 
SNela is considered. In Sect. 2 we introduce the model, in Sect. 
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3 we give our definition of refill timescale, in Sect. 4 we present 
our results and finally in Sect. 5 some conclusions are drawn. 

2. The model 

We simulate a SF region and the late evolution of the hot cavi- 
ties by means of a 2-D hydrodynamical code in cylindrical co- 
ordinates, described in detail in Paper I and references therein. 
The set-up of the model is taken from Paper I, namely it re- 
produces the main features (H i mass and distribution, total dy- 
namical mass) of the dirr galaxy NGC1569. In particular, the 
distribution of gas is approximately ellipsoidal, with a ratio be- 
tween minor and major axis of ~ 0.5, in agreement with ob- 
servations (Reakes 1980). Such a flattened initial distribution 
of gas favours the development of a galactic wind in the polar 
direction, where the pressure gradient is steepest. Although the 
aim of this paper is not to simulate the detailed chemical and 
dynamical evolution of this specific object (as done in Paper 
I), NGC1569 is a good and well-studied example of a post- 
starburst galaxy (Israel 1988), therefore a good benchmark to 
study the refill process after an episode of SF. In order to de- 
velop a physical feeling for the dependence of this process, we 
test two possible values of the total Hi mass: ~ 10^ Mq (mod- 
els labeled "L") and ~ 1.8 x 10** Mq (models labeled "H"). 
According to the works of Angeretti et al. (2005) we concen- 
trate on the SF occurring in the central part of the galaxy, al- 
though we do not intend to reproduce in detail the SF history 
of NGC1569. We consider single SF episodes of different du- 
ration and intensity. We consider either models in which the 
SF lasts 25 Myr at a rate of 0.5 Mq yr"' (second label "S"), or 
models in which the duration of the SF is 200 Myr at a rate of 
0.05 Mq yr ' (second label "L"). Therefore, for instance, the 
model labeled "LS" is characterized by a total H i mass of ~ 
10** Mq and a SF rate of 0.5 Mq yr"' lasting for 25 Myr. Model 
parameters are summarized in Tab. ^ 

The SNela rate is calculated according to the so-called 
Single-Degenerate scenario (namely C-O white dwarfs in bi- 
nary systems that explode after reaching the Chandrasekhar 
mass because of mass transfer from a red giant companion). 
This rate, in the case of short SF episodes, peaks after ~ 10^ yr 
and then decUnes proportionally to ~ t ' (Greggio & Renzini 
1983; Matteucci & Recchi 2001). The chemical evolution of 
our models is also followed, coupling the hydrodynamical sim- 
ulations with chemical yields coming from SNell, SNela and 
intermediate-mass stars. In agreement with Paper I, the sets of 
yields we adopt for massive and intermediate-mass stars are the 
ones calculated by Meynet & Maeder (2002). Details on how to 
trace the chemical enrichment of gas in this kind of simulations 
can be found in Paper I and references therein. 

3. Definition of refill timescale 

According to Paper I, in this work, all the energy and mass 
input occurs in the central 200 x 200 pc^ of the galaxy. This 
assumption comes from the fact that presumably the SF occur- 
ring in the center of NGC1569 comprises the SF of the whole 
galaxy (see Paper I; Greggio et al. 1998). We will therefore de- 
fine as refill timescale the time elapsed from the end of the SF 



phase to the moment at which the physical size of the hot cav- 
ity becomes smaller than the initial 200 x 200 pc^, namely the 
time after which, in spite of the input of energy from SNela, 
a fraction of the central star forming region has temperatures 
below 2 X lO"* K. 

4. Results 

4.1. A typical evolutionary sequence 

As a representative example of the expected evolution of a 
galaxy as a consequence of the energy input of SNela, we 
consider here model HL, namely the model with a large to- 
tal Hi mass and an episode of SF lasting 200 Myr at a rate 
of 0.05 Mq yr"'. Snapshots of the evolution of this model are 
shown in Fig. ^ As one can discern from the first panel, the 
combined energy of SNela and SNell is able to break-out of 
the galaxy at t ~ 160 Myr and a weak galactic wind is formed. 
The superbubble starts funneling through the H i due to the ini- 
tial stratification of the ISM (see Sect. 2) and to the structuring 
of the gas formed by the ongoing dynamical instabilities. The 
energy supply by SNell ends at t ~ 230 Myr (200 Myr is the 
duration of the SF; ~ 30 Myr is the lifetime of a 8 M© star, 
the smallest able to give rise to a SNII in our models). After 
this time interval, SNela still provide enough energy to sustain 
the outflow (panel 2). The funnel begins to shrink at t ~ 300 
Myr and at ~ 340 Myr the outflow has almost completely dis- 
appeared (panel 3). At ~ 400 Myr the cavity has approximately 
the original size of the SF region (panel 4) and from now on 
it recedes further towards the center The last panel of Fig. ^ 
shows the density contours a t ~ 440 Myr and at this time the 
cavity of hot gas (with temperatures of the order of 10^-10^^ 
K) has a size of approximately 100 x 200 pc^; smaller than the 
SF region. According to the definition given in Sect. 3, the refill 
of the cavity has occurred. This result might however depend 
on the dispersion of SNela progenitors (not taken into account 
in this work). Assuming a velocity dispersion of a few km s"' 
the progenitors can travel a few hundreds pc after ~ 100 Myr 

The inflow velocity of the cold gas toward the center of the 
galaxy is modest (of the order of 5-10 km s"') (approximately 
the sound speed of the cold gas), therefore the ram pressure as- 
sociated with this flow is of the order of p^^^ = Pinfaii • ^^ifaii ~ 
10"'^ erg cm"^. This value is three to five orders of magnitudes 
smaller than the ram pressure associated with the inflows sim- 
ulated by Tenorio-Tagle & Munoz-Tunon (1997). The inflows 
considered by these authors should therefore have a different 
origin, such as the infall of a high-velocity cloud directly into 
the center of a galaxy or a violent large-scale disturbance in a 
galactic potential due to a close encounter (Larson 1987). 

It is also worth pointing out that the galactic potential well 
is dominated by a quasi-isothermal dark halo with a core radius 
of 1 kpc (see Paper I for details), therefore the gravitational ac- 
celeration in the center of the galaxy is small. The infall of cold 
gas along the disk is then mainly driven by the pressure gra- 
dient originating from the pressure loss of the cavity once the 
supershell breaks out (MacLow & McCray 1988). However, 
one has to emphasize that in these simulations, at variance with 
similar studies considering only SNell as a source of energy 
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Fig. 1. Density contours and velocity fields for model HL (see Tab.Q at five different epochs (evolutionary times are labelled in 
the box on the upper right corner of each panel). The logarithmic density scale is given in the strip on top of the figure. In order 
to avoid confusion, velocities with values lower than 1/10 of the maximum value (indicated for each panel in the upper right box) 
are not drawn. 



(e.g. DB99), there is no buoyancy of the hot gas at the end of 
the SF, since SNela provide continuously energy. 



plenishment of the SF region close to 1 Gyr in the case of a 
prolongated SF (model "LL"). 



4.2. The evolution of the other models 

In the last column of Tab. ^ we report of the refill timescale 
for each of the considered models. The refill timescale, as ex- 
pected, strongly depends on the total mass of Hi gas at the 
beginning of the simulation. For the extreme model HS (high 
total H I mass, short SF), the refill occurs already ~ 100 Myr af- 
ter the SNell ceased. The tabulated refill timescales associated 
with the "H" models are consistent with the ones derived for 
the standard model studied by DB99, in spite of the fact that no 
SNela were considered in their simulations. This is due to the 
fact that their assumed luminosity in the starburst phase (which 
lasts 30 Myr in their models) is ~ 10 times larger than the one 
adopted in our simulations. DB99 also considered a model with 
a weaker starburst (their model SBl) and in this case the cavity 
begins to shrink immediately after the end of the energy input 
phase and the cold gas reaches the origin at t ~ 70 Myr, much 
earlier than in our simulations. The models characterized by a 
lower H i mass show larger refill timescales, which bring the 
theoretical interval between the beginning of the SF and the re- 



4.3. The impact of a new episode of SF on a refilled 
cavity 

As demonstrated in Recchi et al. (2004) and in Paper I, if a new 
episode of SF occurs in a hot cavity, the metals newly synthe- 
sized in the ongoing starburst are not detectable by the optical 
spectroscopy since they are either directly channelled through 
the galactic funnel, or they are found in a too hot medium. This 
is not the case when the gap between the two episodes of SF is 
large enough to allow the refill of the cavity. In this case, the 
starburst injects metals in a much colder and denser medium, 
where thermal conduction, thermal instabilities and eddies al- 
low a fast cooling of the metals and an effective mixing with 
the surrounding unpolluted gas. In order to show this, we run a 
model in which, starting from the end of model HS (see Tab.Q) 
we add a second burst of SF, with the same characteristics as 
the first one, namely, a SF duration of 25 Myr with a rate of 0.5 
Mq yr"' . The onset of this new SF occurs at ~ 200 Myr, when 
most of the central 200 x 200 pc^ of the galaxy are filled with 
cold gas. 
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Table 1. Model parameters and refill timescales 



Model" 


Hi mass (10'* Mq) 


SF duration (Myr) 


SF rate (Mq yr"') 


Refill timescale (Myr) 


LS 


1. 


25 


0.5 


415 


LL 


1. 


200 


0.05 


600 


HS 


1.8 


25 


0.5 


125 


HL 


1.8 


200 


0.05 


200 



"The identification of a model is made through the notation XY, where X indicates the gas mass ("L" for low and "H" for high) and Y indicates 
the duration of the SF episode ("S" for short and "L" for long). 
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Fig. 2. Evolution of 12 + log (O/H) for model HS (see Tab.[D 
with a second episode of SF at t = 200 Myr. 

In Fig. 12 we plot the evolution of 12 + log (O/H) for this 
model, until an evolutionary time of ~ 240 Myr As one can 
see in this plot, the chemical composition changes significantly 
~ 10-15 Myr after the onset of the last burst of SF. The oxy- 
gen composition increases suddenly as a consequence of the 
input of freshly produced metals by the ongoing SF. The mix- 
ing timescale of the metals produced in this last episode of SF 
is therefore of the order of 10-15 Myr, in agreement with what 
is found by Recchi et al. (2001) in the case of a single instanta- 
neous burst. 

5. Conclusions 

In this paper we have analyzed the late evolution of a galaxy af- 
ter a single SF episode in order to explore the timescale needed 
to refill the cavity in the center of the galaxy with cold gas 
once the SF has ceased. This timescale strongly depends on the 
amount of gas initially present inside the galaxy and ranges 
between 125 (model with a large initial amount of Hi and 
short and intense SF) and 600 Myr (model with a smaller ini- 
tial Hi mass and milder and long SF). This means that large 
Hi holes (like the ones observed in many dirr galaxies) can 
survive a few hundred Myr after the last OB stars have died. 
The refill of the cavity is mostly due to the pressure gradient 
created after the superbubble breaks out the disk. 

A SF occurring in the refilled cavity would produce met- 
als which mix with the surrounding unpolluted medium in a 



timescale of the order of 10-15 Myr, at variance with what 
happens if the center of the galaxy is still occupied by hot and 
diluted gas. 
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